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We review the construction and applications of exactly Poincare invariant quantum mechanical 
models of few-degree of freedom systems. We discuss the construction of dynamical representations 
of the Poincare group on few-particle Hilbert spaces, the relation to quantum field theory, the 
formulation of cluster properties, and practical considerations related to the construction of realistic 
r S^ , interactions and the solution of the dynamical equations. Selected applications illustrate the utility 

Oh of this approach. 
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PACS numbers: 24.10.jv and PACS ll.80.-m 



I. INTRODUCTION 



While there is strong evidence that QCD is the theory of the strong interactions, direct calculations of scattering 
iy-j ' observables in QCD with mathematically controlled errors are difficult at some important energy scales. These 
difficulties are particularly significant at the few-GeV scale, where perturbative methods are not applicable. This is 
C\| ■ an interesting energy scale because it is the scale where sensitivity to sub-nuclear degrees of freedom is expected to 
I/"*) | begin. Mathematical models that are motivated by QCD may provide useful insight into the dynamics at these energy 
. scales. 

Poincare invariant quantum mechanics is one of a number of approaches that can be used to model systems of 
. strongly interacting particles at the few GeV energy scale. At the simplest level it is quantum mechanics with an 
underlying Poincare symmetry. While Poincare invariant quantum mechanics can be treated as a phenomenology 
, that is independent of QCD, it can also be related to QCD. Poincare invariant quantum mechanics has proved to be 
useful in applications, but there are no textbook treatments of the subject. 

Historically, Poincare invariant quantum mechanics was first articulated by Wigner who pointed out that a 
necessary and sufficient condition for a quantum theory to be relativistically invariant is the existence of a unitary 
ray representation of the Poincare group on the quantum mechanical Hilbert space. Wigner's work did not have a 
significant impact on applications of quantum field theory, but it directly motivated attempts to provide an axiomatic 
[2hJ] foundation for quantum field theory. These axioms provide a Hilbert space formulation of quantum field theory 
that can be directly related to Poincare invariant quantum mechanics. 

Dirac [|[ studied the problem of constructing the Poincare Lie algebra for systems of interacting particles. He 
observed that the presence of interactions in the Hamiltonian implied that at most a sub-algebra of the Poincare Lie 
algebra could be free of interactions. He identified the three largest sub-algebras, and classified dynamical models 
according to which sub-algebra remained free of interactions. Dirac used the terms instant, point, and front-forms of 
dynamics to label the different kinematic sub-algebras. Bakamjian and Thomas [|| provided the first construction of 
the full Poincare Lie algebra for a system of two interacting particles in Dirac's instant-form of the dynamics. Coestcr 
Q generalized Bakamjian and Thomas' construction to systems of three interacting particles. His construction also led 
to a S matrix that satisfied spacelike cluster properties. Sokolov Q provided a complete construction of the Poincare 
Lie Algebra for a system of N interacting particles in Dirac's point-form of the dynamics that was consistent with a 
stronger form of spacelike cluster properties, where the Poincare generators satisfy cluster properties. This stronger 
form of cluster properties provides a simple relation between the few and many-body systems that is difficult to realize 
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in theories satisfying only S-matrix cluster properties. Coester and Polyzou Q provided the complete solution for 
systems of TV-particles in all three of Dirac's forms of the dynamics satisfying the strong form of cluster properties. 
A more general construction based on only group representations, that has Dirac's form of dynamics as special cases, 
was given in [13, dl- The subject was reviewed by Keister and Polyzou in [l2T |. 

There have been many applications of Poincare invariant quantum mechanics in all three of Dirac's forms of 
dynamics. The earliest applications involved the study of electromagnetic probes on mesons, nucleons, and nuclei. 
Some of the relevant papers are [1314271. T he first three-nucleon bound state calculation using this framework was 
performed by Glockle, Coester and Lee [28j . Calculations of the triton binding energy with realistic interactions have 
been performed recently (29j . Applications to nuclear reactions appear in (30l - [3a | which include reactions with particle 
production [36| . 

This mini-review is limited to theories that are formulated by constructing exact unitary representations of the 
Poincare group on few-particle Hilbert spaces. There are many other approaches to relativistic quantum mechanics 
that have been successfully applied at the few GeV scale. Each one emphasizes different desirable features of the full 
field theory, however when the number of degrees of freedom is limited, it is impossible to satisfy all of the axioms of the 
underlying field theory. Our preference for using Poincare invariant quantum mechanics is based on three observations: 
(1) many computational methods successfully used in non-relativistic quantum mechanics can be directly applied in 
Poincare invariant quantum mechanics, (2) the theories involve a finite number of degrees of freedom, allowing exact 
numerical calculations of model predictions, (3) the theories share most of the axiomatic properties of quantum field 
theory and there is a direct relation to the Hilbert space formulation of field theory. The fundamental property of the 
quantum field theory that is given up in order to have a theory of a finite number of degrees of freedom is microscopic 
locality. The justification for this choice is that microscopic locality is not an experimentally testable property since 
probing a system at arbitrarily short distance scales requires arbitrarily large energy transfers. In addition, Poincare 
invariant quantum mechanics does not have a large enough algebra of observables to localize particles in arbitrarily 



small spacctime regions. One manifestation of this is the absence of a reasonable position operator [3JJ in relativistic 
quantum theories of a finite number of degrees of freedom. 

In the next section we discuss the construction of representations of single-particle Hilbert spaces. In section ?? we 
discuss irreducible representations of the Poincare group that act on the single-particle Hilbert spaces. In section ?? 
we construct a dynamical representation of the Poincare group by adding interactions to the mass Casimir operator of 
a non-interacting irreducible representation constructed from tensor products of single particle representations. The 
strong and weak form of cluster properties are discussed in section ??. The formulation of the three-body problem 
is discussed in section ??. The relation to quantum field theory is discussed in section ??. Selected few-nucleon 
applications arc discussed in section ??. 



II. PARTICLES, HILBERT SPACES AND IRREDUCIBLE REPRESENTATIONS 

Experiments measure observables that describe the state of free particles by considering how the particles interact 
with classical electromagnetic fields. A complete experiment measures the linear momentum and spin state of each 
initial and final particle. There is a natural connection with these single-particle observables and irreducible repre- 
sentations of the Poincare group. The Poincare group has ten infinitesimal generators. These Hcrmitian operators 
include the Hamiltonian which generates time translations, the linear momentum operators which generate space 
translations, the angular momentum operators which generate rotations, and the rotationless boost generators which 
generate transformations that change the momentum of the particle. From these ten infinitesimal generators it is 
possible to construct two Casimir invariants, four independent commuting Hcrmitian observables and four conjugate 
operators. The Casimir invariants fix the mass and spin of the particle. Eigenvalues of the commuting observables 
label the states of the particle, and the conjugate operators determine the spectrum of the commuting observables 
and thus the allowed states of the particle. 

For a standard description of a particle, the commuting observables can be taken to be the three components of the 
linear momentum, and a component of a spin operator. The spectrum of the momentum is R 3 , while the spectrum of 
a component of the spin vector takes on discrete values in integer steps from —j to j. In this case the Hilbert space is 

Hmj =L 2 (R 3 )®C 2j+1 . (2.1) 

Single-particle states are represented by wave functions, ip(jp,fi) = ((to, j)p, fJ,\ip). The relations of the operators 
to, j 2 , p and j • z to the Poincare Lie Algebra determines a unitary representation, U\ (A, a), of Poincare group on H m j ■ 



((77l,j)p,/i|l7l(A,o)|V') = 
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where the Poincare group Wigner function is 



S(p Ap')J ^ge^ [R wc (A,p)] (2.3) 
y v m (p ) 

and R wc (A,p) is a Wigner rotation. 

Because a sequence of Lorentz boosts that start and end at the rest frame generally define a rotation, in order 
to obtain an unambiguous definition of a spin vector for all values of the particle's momentum, it is necessary to 
define a standard way to measure a spin observable. The above representation implicitly defines the spin projection 
by its value in the particle's rest frame after the particle is transformed to the rest frame with a rotationless Lorentz 
transformation. This is one of an infinite number of possible choices of spin observables. This choice is consistent 
with the "canonical" spin that appears in standard Dirac u and v spinors. Different spin observables are related 
by momentum-dependent rotations that lead to different couplings to the electromagnetic field. This ensures that 
measurable physical quantities are independent of the observables used to label single particle states. 

These single-particle representations are irreducible, and all positive-mass positive-energy irreducible representa- 
tions of the Poincare group can be put in this general form. These irreducible representations will be important in 
formulating dynamical models. In general, any unitary representation of the Poincare group can be decomposed into 
a direct sum or direct integral (for continuous mass eigenvalues) of irreducible representations. We will build the 
dynamical unitary representation of the Poincare group out of the non-interacting irreducible representations. 

III. POINCARE GROUP WIGNER FUNCTIONS AND KINEMATIC SUBGROUPS 

In the previous section we represented single-particle wave functions in the basis of generalized eigenstates 
|(m,j),P,M> • 

The state of the particle could be also alternatively determined by measuring the particles' four velocity, = 
(y/l — v • v, v), and spin projection: 

\{m,j),w,n) = \(m, j),p(v,m),^)m 3/2 , (3.1) 

or the light-front components of the four momentum p + = y m 2 + P 2 + P ■ z, p± = (p ■ x, p ■ y), and light-front spin 
projection: 

\(m,j),p + ,p ± ,») = £ \(rn,j)Mp + ,P±,rn)^')J^^D^^B- 1 (p)B f (p)} (3.2) 

A*'=— J ' 

where B~ 1 (p)Bf(p) is a Melosh rotation [38| . defined by a light-front-preserving boost followed by the inverse of a 
rotationless boost. The different basis choices are related to the basis |(m,j),p,/x) by the unitary transformations in 
(|3.ip and p. 21) . The light-front preserving boosts have the desirable property that they form a group, which means 
that there are no Winger rotations for any sequence of light-front preserving boosts that start and end in the rest 
frame; the price paid for this is that the Wigner rotation of a pure rotation is not equal to the rotation. 
The Poincare group Wigner functions depend on the choice of basis. The Wigner functions 

Ki^A^ ~ (Kj)v,A^i(A,a)|(m,j)v', A (3.3) 

^.^.p-ylA,"] := ((m,j)p + ,P±,M|C/i(A,a)|(m, J y+,p' ± , / i') (3.4) 

are related to the Wigner function (|2.3p by the unitary transformations (|3.1|) and (|3.2[) . 

While the concept of a kinematic subgroup does not make sense for a single particle, the kinematic subgroup for 
an instant-form dynamics is the subgroup of the Poincare group that leaves the Wigner function (|2.3[) independent of 
mass; the kinematic subgroup for a point-form dynamics is the subgroup of the Poincare group that leaves the Wigner 
function (|3.3[) independent of mass; the kinematic subgroup for a front-form dynamics is the subgroup of the Poincare 
group that leaves the Wigner function (|3.4p independent of mass. Different mass-independent subgroups appear in 
different irreducible bases because the unitary transformations relating the irreducible bases (|3 . 1 [) and p.2[) to the 
basis \(m,j)p,n) depend on the particles' mass. These mass-independent subgroups become kinematic subgroups in 
dynamical models because the mass acquires an interaction while the other operators used to construct dynamical 
irreducible bases remain interaction free. More generally, it is possible to define perfectly good single-particle bases 
where the identity is the only subgroup where the corresponding Wigner function is independent of mass. 
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IV. TWO-BODY MODELS - CLEBSCH-GORDAN COEFFICIENTS 

The two-body Hilbert space is a tensor product of two singlc-particlc Hubert spaces, H = H mi j 1 <8> 7~Lm 2 j 2 ■ The 
non-interacting representation of the Poincare group on H is the tensor product of two single-particle (irreducible) 
representations of the Poincare group. C/ (A,a) := Ui(A,a) ® U 2 {A,a) . While the single-particle representations of 
the Poincare group are irreducible, their tensor product is reducible. Formally the tensor product representation can 
be expressed as a direct integral of irreducible representations, 

U {A,a) = y2 @ dmU , mtj ,i t3 (A,a), (4.1) 



E 



where Uo, m ,j,i,s(A,a) are mass m spin j irreducible representations of the Poincare group. The quantum numbers / 
and s are invariant degeneracy parameters that distinguish multiple copies of the irreducible representations with the 
same m and j. They have the same quantum numbers as the spin and orbital angular momentum. The mass m is 
the two-particle invariant mass that has a continuous spectrum starting from mi + m 2 . The Poincare group Clcbsch- 
Gordan coefficients relate the tensor product representation to the direct integral of irreducible representations and 
satisfy 

J V p,i: P >,n> I A ' °]dp'{{m,j, I, s), p', //|(mi, ji), pi, Hi; (m 2 ,j 2 ), P2, H2) = 

E / '(( m > S )' P ' '"K"' 11 ' Jl)iP'li Ml 5 ( TO 2,j2),P 2 ,M2> rf PlP2 X 

V mu3 l \A,aW m , 2 ^ [A, a]. (4.2) 

The Clebsch-Gordan coefficients, ((m, 7, d), p, ft|("T-i v 7i ) , Pi,fii; (?7i2, J2), P2, ^2), d := {£, s}, are basis-dependent and 
are known in all three of the representations ( [3, LL2], [30, |40( ) . 

The two-body irreducible basis states look similar to relative and center of mass variables in non-relativistic quan- 
tum mechanics; but they differ in the structure of the Poincare group Clebsch-Gordan coefficients, which contain 
momentum-dependent spin rotation functions and non-trivial kinematic factors that ensure unitarity. 

The basis states {\(m,j,l,s),p,fi)} transform irreducibly under {/ (A, a): 

U (A,a)\(m,j,d),p,»)=J2 j dp'\(m,jJ),p\»')V™^,. p JA,a}. (4.3) 

While (|4.3[) is not the dynamical representation of the Poincare group, by working in this non-interacting irreducible 
basis it is possible to construct dynamical representations by adding an interaction v, which in this basis has a kernel 
of the form 

((m'jVO.P'.M'MKi.^P.M) = <Sj' J ^ M *(p , -p)(m',d'||« J *||m,d), (4.4) 

to the non-interacting two-body mass operator . This in tera ction has the same form as a typical Galilean invariant 
non-relativistic interaction if we replace m = y m \ + k 2 + m\ + k 2 by k and d = {I, s} by (I, s). 

We define the dynamical mass operator M := \J m\ + k 2 + \J m| + k 2 + v . Simultaneous eigenstates of M, p, j 2 
and j • z can be constructed by diagonalizing M in the irreducible non-interacting basis. These eigenfunctions have 
the form 

((k',f, I', a'), p', M'|(A,i), p, n) = fyj Vm*(p' - P)^d ( k ^> s ) ( 4 - 5 ) 
where the wave function, <fi\.j(k 2 , l,s), is the solution of the mass eigenvalue problem with eigenvalue A: 



(A - J m\ + k 2 - ^Jrn\+k 2 )(t>\ j (k, I, s) 



/ k' 2 dk'^2 E (*,i,*|V J |*',I / I s , )0A3(A',i' J a')- ( 4 - 6 ) 

J ° s' l'=\js\ 
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The dynamical unitary representation of the Poincare group is defined on this complete set of states, |(A,_j), p, n) 

by 

[/(A,a)|(A,j),P,M>= £ jdp'KXJlp'^'^^JKa}. (4.7) 

The relevant dynamical feature is that the Poincare group Wigner function now depends on the eigenvalue A of the 
dynamical mass operator, which requires solving (|4.6[) . 

Because of the choice of basis, the Poincare group Wigner function P p ; J ' , [A, a] has the same structure as the 
Wigncr function (|2.3|) and thus has the property that when (A, a) is in the three-dimensional Euclidean subgroup, 
it is independent of the mass eigenvalue A, which means that for this dynamical model the kinematic subgroup is 
dictated by the choice of representation used to define the irreducible basis. 

Even though the dynamics has a non-trivial interaction dependence, it is only necessary to solve (|4.6[) . which is 
analogous to solving the center of mass Schrodingcr equation in the non-relativistic case. 

This construction can be repeated using different irreducible bases, such as (|3 . 1 [) or (|3.2p . where the Wigner functions 
have different mass-independent symmetry groups. For these bases if we choose to use the interactions 

((k',j',l',s'),v',ti'\v point \(k,j,l,s),v,iJ,) = S rj S^^S(V -v)(k',l',s'\\v j \\k,l,s) (4.8) 
((k',f, I', s'),p' + , p' ± ,(i'\v fron t\(k,j, I, s),p + ,p ± ,fj) = 

fyiV^(Pi - P±)S(P' + - P + )(k', I', s'lkUfc, I, s), (4.9) 

where the reduced kernels (k', I', s'\\v J \\k, /, s) arc the same in (|4~4)) . (|4~8)l , and (|4~9)) in the bases (|2~2|) . (j3~T]) . and (|3~2l) . 
respectively, and construct dynamical eigenstates of the form 

|(A,j),v,/i), \(X,j),p + , P± ,fi) , (4.10) 

then equation (|4.6[) still determines the binding energy and scattering phase shifts. It follows that the resulting 
two-body models have the same bound-state and scattering observables, however each of the resulting unitary repre- 
sentations of the Poincare group has a different kinematic subgroup. The dynamical irreducible eigenstates transform 
like |(A,_7'), p, fi) with the Wigncr function (|4.7[) replaced by (|3.3j) or (|3.4p where m is replaced by A. This makes these 
unitary transformations dynamical. 

The mass operators and interactions, v, v point and Vf ront are distinct operators, but the three representations 
are related by unitary transformation that leave the binding energies and scattering observables unchanged. The 
dynamical calculations are identical in all three cases and are given by solving (|4.6[) . This shows that dynamical 
models with different kinematic subgroups are equivalent and cannot be distinguished on the basis of any experimental 
observations. 



V. CLUSTER PROPERTIES - EKSTEIN'S THEOREM 

An important feature of non-relativistic quantum mechanics is that the same interactions appear in the few and 
many-body problems. Specifically, the Hamiltonian becomes a sum of subsystem Hamiltonians when the short- 
ranged interactions between particles in different subsystems are turned off. In the relativistic case the corresponding 
requirement is that the unitary time-translation group breaks up into a tensor product of subsystem groups when the 
system is asymptotically separated into independent subsystems. We call this the strong form of cluster properties. 

The observable requirement is that the S'-matrix clusters. We call this the weak form of cluster properties because 
it follows from the strong form of cluster properties, however because different Hamiltonians can have the same 
S'-matrix, the weak form of cluster properties does not imply that the same interactions appear in the few and many- 
body Hamiltonians. Because of this, in order to maintain a simple relation between the few and many-body problem, 
we require that Poincare invariant quantum theories satisfy the strong form of cluster properties. 

A theorem of Ekstein [4lj provides necessary and sufficient conditions for two short-ranged interactions to give the 
same S matrix. The requirement is that the Hamiltonians are related by a unitary transformation A satisfying the 
asymptotic condition 

lim \\(I-A)Uo(tm\\=0 (5.1) 
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where Uo(t) is the non-interacting time translation operator. We refer to unitary transformations with this property as 
scattering equivalences. It is important that this condition be satisfied for both time limits; to appreciate the relevance 
of this condition consider two Hamiltonians with different repulsive potentials. Because these Hamiltonians have the 
same spectrum and multiplicities they are related by a unitary transformation, however the derived S-matrices may 
have different phase shifts. The phase shifts differ if and only if two time limits do not agree. 

Scattering equivalences that preserve weak cluster properties but not strong cluster properties exist and are the key 
to restoring the strong form of cluster properties in Poincare invariant quantum theory. The strategy is illustrated in 
the formulation of the three-body problem in the next section. 

VI. THREE-BODY PROBLEM 

The strong form of cluster properties implies that given a set of dynamical two-body generators, the three-body 
generators necessarily can be expressed as sums of one, two and three-body operators 



H — Hi + H-2 + Hs + H\2 + H-23 + H31 + H123 (6-1) 

P = Pi + P 2 + P 3 + P12 + P23 + P31 + P123 (6.2) 

J = J : + J 2 + J 3 + J 12 + J 23 + J31 + J123 (6.3) 

K = Ki + K 2 + K 3 + K12 + K 23 + K31 + K 123 . (6.4) 



The one and two-body operators in (|6.H6.4p are the same operators that appear in the two-body problems, while the 
three-body operators, i?i 23 , Pi 23 , Ji 23 and Ki 23 are the only new ingredients in the three-particle generators. 

It is easy to show that if the generators have this form it is impossible to satisfy the Poincare commutation relations 
if all of the three-body operators vanish. However, although the commutation relations put non-linear constraints on 
these operators, it will become clear that the solutions are not unique. 

To avoid solving the non-linear problem of satisfying the commutation relations, it is more productive to start by 
first satisfying the commutation relations at the expense of strong cluster properties. This can be done by applying 
the method of section ?? directly to the three-body problem. This involves adding suitable interactions to the 
non-interacting invariant three-body mass operator. 

To begin the construction we consider a three-body system where only one pair of particles interact. The relevant 
basis is a non-interacting three-body irreducible representation of the Poincare group. It is constructed by successive 
pairwise coupling using the Poincare group Clebsch-Gordan coefficients. If we assume that particles one and two are 
the interacting pair then preferred order of coupling would be (12) — > ((12)(3)): 

|pi, Mi) ® |P2,M2) \(ki2,li2,si2, iia)pi2, M12) (6.5) 
\(kl2, I12, S12, ji 2 )pi 2 , M12) ® |p 3 , 1^3) -> 

\(Q, £(12)(3),<S(12)( 3 )./(12)(3), fcl2, Zl2, »12, jl2)P, fi) ■ (6-6) 

We introduce the following shorthand notation for the basis states in these equations. We write (|6.5[) as |1(8>2) — > | (12)) 
and (16.61) as |(12)®3) — s- |(12)(3)). Using this notation we define two different embeddings of the two-body interaction 
in the three-body Hilbert space using the two representation in (|6 .€>[) : 

((12)'®3>i2®|(12)®3) = 

( fc 12^12 I s i2lk J |l fc 12^12,Sl 2 ,)J(pi 2 -Pi 2 )(5(P3 ~ P3)Sj' 12 j l2 Sj' 3 j 3 (6.7) 

and 

((12)'(3)> 12 |(12)(3)> = 
(fci 2 ,/ 12 ,s' 12 ||w- 7 ||fci 2 ,Zi 2 ,si 2 ,)(5(p' - p)x 
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8(<f-q) x x X x x (p,r\ 

j;i2)(3)J<12)(3) ^2il2 i' (12 , ( 3) i (12)(3) S; i2)(3) S (12)(3) 0MV 



<l 2 



where the reduced kernel, (k' 12 , l\ 2 , s 'i2\\ v:> \\k12, Z12, Si2> ), is identical in (|6.7[) and (|6.8|) . These expressions define 
different interactions, («i 2 ® 7^ 1*12) • 

We use these two interactions to define two different 2 + 1-body mass operators Af^)®^) and M( 12 )( 3 ) defined by 



\/(\A\/ TO i + + yjml + k\ 2 + Wl20 ) 2 + p^) + ^jw? +^§) 2 - p 2 ) (6.9) 




(6.10) 



Because of the invariance principle [42h44| the S'-matrix can be computed by replacing the Hamiltonian by the mass 
operator (this is equivalent to evaluating the ^-matrix in the three-body rest frame) in the standard time-dependent 
representation of the scattering operator. 

M( 12 )g,(3) is the mass operator of the tensor product of a two-body representation involving particles one and 
two and a spectator representation of the Poincare group associated with particle three, Ui 2 (A,a) ® t/3 (A, a). By 
construction it is consistent with the strong form of cluster properties. The mass operator Af( 12 )(3) commutes with the 
three-body spin and commutes with and is independent of the total three-body momentum and 2-component of the 
three-body canonical spin. Simultaneous eigenstates of M( 12 )(3), p, j 2 , j z are complete and transform irreducible with 
respect to the Poincare group. This defines a dynamical unitary representation of the Poincare group, C/( 12 )(3)(A, a), 
on the three-body Hilbert space following the construction of section ??. 

The scattering operators associated with both of these operators are related by 

((12)®(3)|S {12)fi9(3) |(12)®(3)> = 



(KzJli^lsWS 3 ||&12,Jl2, Sl2,)6(p' 12 -P12WP3 - P3)Sj' 12 j 12 Sj' 3 j 3 (6.11) 

and 



<(12)(3)|S (12)(3) |(12)(3)) = (k' 12 , l' 12 , s' 12 \\S^\\k 12 ,l 12 ,s 12 ,)6{p' - p)x 



q 2 ^';i2)(3)J(12)<3)^';2J12^' ( i 2 ) ( 3) i <12)(3)^S[ 12)(3) S ( i 2 ) ( 3)^'A' (6-12) 

where the reduced two-body kernels (k[ 2 , l'i 2 , s' 12 \\S 3 \\ki 2 , l\ 2 , s\ 2 , ) are identical. Because the delta functions become 
equivalent when they are evaluated on shell, the S matrices in both representations are identical. Ekstein's theorem 
implies the scattering equivalence j4(i2)(3)^12(A, a) <g> Ua(A, a)j4| 12 w 3 , = f/(i2)(3)(A, a). 

To construct a dynamical representation of the Poincare group with all three particles interacting we first construct 
the mass operator 

M = Af(ia)( 3 ) + M( 23 )(i) + M (31) - 2M . (6.13) 

Because each term in (|6 . 1 3|) commutes with p,j 2 ,j z , and is independent of p and j z it follows that simultaneous 
eigenstates of M,p,j,j z are complete and transform irreducibly, thus defining a dynamical unitary representation, 
U (A, a), of the Poincare group on the three- nucleon Hilbert space. 

Because each of the 2 + 1 mass operators in (|6.13[) is scattering equivalent to 2 + 1 mass operators associated with 
a tensor product representation, it follows that M can be expressed as 



M = A(12)(3)Af(i2)®(3)A (12)(3) + A(23)(l)Af(23)®(l)^ (23)(1 



^(31)(2)M(31)®(a)4»l)(3) " 2M °- (6 ' 14) 

From this representation it follows that when interaction between the i th particle and the other two particles are 
turned off that 



U (A, a) -> A m{i) U m (A, a) ® U l {K, a)A ] (jm) 



(6.15) 



8 



which formally violates the strong form of cluster properties. 

The strong form of cluster properties can be restored by transforming U(A,a) with the product = 
^(i2)(3)^(3i)(2)^(23)(i)- Because products of scattering equivalences are scattering equivalences, this does not change 
the three-body S matrix. This transformation also restores strong-cluster properties, because — > A^^,^ when 
the interactions between particle i and the other two particles are turned off, canceling off the extra unitary trans- 
formations in (|6.15[) . The undesirable feature of A is that the individual A^^^s do not commute, so it introduces 
an exchange asymmetry that does not affect the S'-matrix. The exchange symmetry can be manifestly restored by 
replacing the product of the Aq^^'s by a symmetrized product, such as: 

A := e 1 "^ 12 '* 3 '^ 1 "^ 23 )' 1 )^ 1 "^ 31 " 2 )) (6.16) 

E/®(A, a) = A^U(A, a) A. (6.17) 

Equation (|6.17|) defines a unitary representation, C/®(A, a), of the Poincare group that satisfies the strong form of 
cluster properties because 

A A (jmi) (6.18) 
when the interactions between particle i and the other two particles are turned off. Thus 

tf»(A,a) -> A\ jk)(i) U (jkm {K,a)A (jk)(l) = 

A \ 3 km A m(i) u m(^ a ) ® ^(^"K'nt^ww = 

U(j k )(A,a) ® Ui(A,a) (6.19) 

This property ensures that the infinitesimal generators have the additive form (|6.HI6.4j) and (|6.17[) generates the 
required three-body interactions. 

Because there are many other ways to construct symmetric products of non-commuting operators and because it 
is possible to add a three-body interaction to M that commutes with and is independent of the total momentum 
and spin, it is clear the three-body parts of the generators that are required to restore the commutation relations 
are not unique. It is also important to note that it is not possible to use the freedom to add three-body interactions 
to eliminate the three-body interactions required to restore the commutation relations; in this representation the 
generated three-body interactions do not commute with the non-interacting spin. 

This construction can be extended to formulate dynamical models satisfying the strong form of cluster properties 
for any fixed number of particles, isobar models in any of Dirac's form of dynamics. It is even possible to treat 
production beyond isobar types of models. 

Models with different kinematic subgroups can be constructed by starting with different irreducible bases (|3.1[) . 
(|3.2j) . As long as the reduced kernels of the interactions are identical, all of the Bakamjian-Thomas three-body mass 
operators, M , will give identical bound-state and scattering observables. They are related by scattering equivalences 
constructed by applying the unitary transformations 

|(A,j),v,/v.) = |(A,j),p(v,A),/v)A 3/2 , (6.20) 



\(X,j),p+,p ± ,fx, ■■■) = \(X,j), p(p+,p±, A),//, • • .J^-Dl^lB-^Bf^)] (6.21) 

on each invariant subspacc of the associated mass operator. Each of these representation is in turn scattering equivalent 
to a representation that satisfies strong cluster properties and has the same kinematic subgroup. 

Because A is a scattering equivalence, it is only necessary to solve the Faddeev equations for the mass operator M 
in a non-interacting irreducible basis. Furthermore, since all bound state and scattering observables can be computed 
using only the internal mass operator, with the delta functions in p and /i removed, this equation is the same in all 
of Dirac's forms of dynamics when expressed in terms of the kinematic mass and kinematically invariant degeneracy 
quantum numbers. The operators A and the choice of kinematic subgroup are only needed if the three-body system 
is embedded in the four-body Hilbert space or if the eigenstates are used to construct clectroweak current matrix 
elements. 



9 



VII. CONNECTION WITH QUANTUM FIELD THEORY 



Poincare invariant quantum mechanics as formulated by Bakamjian and Thomas resembles non-relativistic quantum 
mechanics more than quantum held theory. The Hilbcrt spaces have the same structure as non-relativistic Hilbcrt 
spaces, the theory is not manifestly covariant, spin 1/2 particles are treated using two-component spinors. In spite of 
these apparent differences there is a direct connection to quantum field theory which we outline below. 

To develop the connection we assume the existence of an underlying quantum field theory with a Poincare invariant 
vacuum and a collection of Hciscnbcrg fields, 4> i {x) 1 where the bold face indicates a multi-component field. The index 
i distinguishes different types of fields. 

In quantum field theory Hilbert-space vectors are constructed by applying functions of smeared fields, 

Mf)= J d 4 xi(x) ■ ^(x) (7.1) 

to the physical vacuum |0). 

Polynomials in the smeared fields applied to the physical vacuum generate a dense set of vectors. The field theoretic 
unitary representation of the Poincare group U^(A, a) acts covariantly on the smeared fields: 

C/ t (A,a)0 4 (f)C/(A,a) = f d i xt(Ax + a)S(A)cf> i (x) (7.2) 



where S(A) is the finite dimensional representation of the Lorentz group appropriate to the field. The covariance of 
the fields implies Poincare transformation properties of test functions that leave the scalar product invariant. 

If the field theory has one-particle states, then there are functions, A, of smeared fields with the property that j4|0) 
is a one particle state. One-particle eigenstates that transform irreducibly with respect to the Poincare group can be 
constructed by projecting A\0) on states of sharp linear momentum and canonical spin. This can be done using the 
unitary representation (|7.2[) of the Poincare group 

\(m,j)p,fj) = A(p,/i)|0) := 

j r 

Y, dRdp°d i xe ipx U{R,x)A\Q)D^{R)5{p 2 +m 2 )6»(P°) (7.3) 

where R is a rotation, dR is the SU (2) Haar measure, U (R, x) is the unitary representation of the Poincare group 
restricted to rotations and spacetime translations, and D^*{R) is a SU{2) Wigner function. 
The normalization of these states can be chosen so 

((m',i')p'>A*'l( m >i)P)M) =tf(p'-p)tfm'm<5j'i<W ( 7 - 4 ) 

It follows from the definitions and the group representation properties that these states transform as mass m spin j 
irreducible representations of the Poincare group: 



j 



U(A,a)\(m,j)p,») = J2 l(™,j)A P ,, 1 ')A»^;[B- 1 (A(p))A5(p)] 1 



< u m (Ap) 

Wm(p) 



£ J dp'\(mJ)p\»')V>™^-, v J^a} 



(7.5) 



To construct scattering states define C(p, fi) := (y 7 m 2 + p 2 A(p, fj,)) — [H, A(p 7 Scattering states are then given 

by the Haag-Ruelle method: dHH] 



|(pi,Mi,--- ,Pn^n) ± ) = Hm C/(-t)n[Q(P J5 ^K lt ^ (Pj) ]|0) (7.6) 
\vm(tc the limits are strong limits after smearing over suitable momentum wave packets. 
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The operators n 7 -[Cj (Pj> Mj)|0) can be considered as mappings from an N particle channel Hilbcrt space, H a , to 
the Hilbert space of the field theory. Vectors in the TV-particle channel Hilbert space are square integrable functions 
in the variables pi,/xi, • ■ ■ , pn, ^n- We denote these operators by ft a ± where a indicates the channel. 

The direct sum of all of the channel Hilbert spaces, including the one-particle channels, defines an asymptotic 
Hilbert space. We define fi± that maps the asymptotic Hilbert space to the physical Hilbert space by 

= ^0 Qi± |f Q! ). (7.7) 

a 

By construction these wave operators satisfy the intertwining relations [46j | 

U(A,a){l± = n±® a U a (A,a). (7.8) 

The Poincare invariant S operator of the field theory is given by 

S = n{ 0_ (7.9) 

where each U a {A,a) is a tensor product of single particle irreducible representations of the Poincare group on the 
channel subspace H a 

Poincare invariant quantum mechanics formulated in the previous sections has the same basic structure. The 
primary difference is that the asymptotic Hilbcrt space for the field theory has an infinite number of channels and 
describes physics at all energy scales, while the Poincare invariant quantum mechanical wave operators involve only 
a subset of these channels that are experimentally relevant only up to a given energy scale. 

If n is a Poincare invariant projection operator on the asymptotic subspaces corresponding channels of a Poincare 
invariant quantum model that also limits the maximum invariant mass of the asymptotic states, then the following 
operator 

w = 0/+no^ m+ = si f _iml m _ (7.10) 

maps an invariant subspace of the quantum mechanical Hilbert space to an invariant subspace of the field theory 
Hilbcrt space in a manner that satisfies 

ns* 9m n = us f u (7.11) 



WU P (A, a) = U f (A, a)W . (7.12) 

These mappings define the relevant relation between the Poincare invariant quantum theory and the underlying field 
theory. 

Thus, for asymptotic scattering states in the range of n the Poincare invariant quantum mechanical theory can be 
designed to give identical results to the field theory. Obviously the two theories differ on asymptotic states that are 
not in the range of H 

Even though the Poincare invariant quantum theory does not satisfy microscopic locality, we see that it can give 
the same S matrix elements as the full field theory at a given energy scale. 



VIII. FEW NUCLEON APPLICATIONS 



In this section we discuss an illustrative set of applications to few nucleon problems. A realistic nuclcon-nucleon 
interaction is needed for these applications. The invariant mass operator for two free nucleons can be expressed in 
terms of a relative momentum as 



=: \J k 2 +m\ + yk 2 +m%. (8.1) 
It is always possible to express the two-body interaction as an addition to k 2 : 



M12 = m i2 + v 12 := J k 2 + 2/j,v nn + m\ + J k 2 + 2fxv nn + ml (8.2) 
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where following [47[ v nn is a realistic nucleon-nuclcon interaction [48J, |49J and \i is the reduced mass 

H = 8.3 

mi + m 2 

In this representation the dynamical two-body mass operator becomes a function of the non-relativistic center of mass 
Hamiltonian: 



M12 = J 2fih nr + m\ + J 2fih nr + m\ (8.4) 

where 

k 2 

h nr = — + V„r- (8.5) 



It is a consequence of the Kato-Birman |42|-|44| invariance principle that the relativistic wave operators for (|8.2[) and 
non-relativistic wave operators for (|8.5[) are identical 



Q nr ± := lim e iHnrt e -iH nr0 t = r lm e ih nr t e -ihot = Hm e iMt e -iM i 
t— >±oo i— >±oo >-±oo 



lim e lM { e- lM o« = i im e ^* e --",ot = i im g^rtg-t^t = ^ ± ( 8 6 ) 

t— >±oo t— >±oo t— »±oc 



where M = M12 in (|8.6[) . The identity (|8.6[) ensures that both scattering operators are identical as functions of k 2 : 

5Vir ^nr+^ m '~ f^.j.^r— S T (8.T) 

Here the relativistic and non-relativistic S are related because the interactions are fit to the same two-body data 
correctly transformed to the center of momentum frame. The non-relativistic Hamiltonian (|8.5[) is NOT the non- 
relativistic limit of (|8.2[) . 

This construction, which first appeared in (47| , shows that existing realistic interactions can be directly used in the 
formulation of a Poincare invariant two-body problem. Equation (|8.4[) implies that the wave functions of (|8.5p and 
(j8.2[) are identical functions of k 2 , 1, s. 

If we replace the interaction in (|4.4[) by the interaction V12 in (|8.2[) and use this in the three-body calculation 
discussed in section ?? then the three-body S'-matrix can be expressed in terms of three-body mass operators: 

Sac = 5 ac - 2mS(M a - M c )T ac (z c ) (8.8) 

which are functions of the transition operators 

T ac (z) = T ac (z) = V C + V a R(z)V c (8.9) 

where a,b,cE {(12)(3), (23)(1), (31)(2)}, 



M {ij)(k) = J (yjm? + k 2 + 2fiv nn + sjm) + k 2 + 2^~ n ) 2 + q 2 + yj ' m 2 k + q 2 . (8.10) 

V a = M a -M V c = J2 Va R{z) = {z-M)~ x (8.11) 

R c (z) = (z- M - V,)- 1 (8.12) 

R(z) = R c (z) + R c (z)V c R(z) (8.13) 

and v nn is the nucleon-nucleon interaction that appears in (|8.2p embedded in the three-nucleon Hilbert space with 
the delta functions in (|6.8|) . The Faddeev equations can be derived using standard methods 

T ab (z) = V b + J2 V cRc(z)T cb (z). (8.14) 
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While it does not make any sense to study the non-relativistic limit of interactions that are constructed by fitting to 
two-body bound and scattering data, we can compare the relativistic and non-relativistic three-body calculations that 
use the same two-body interaction, v nn , as input. In the Poincare invariant quantum mechanics case the Faddccv 
equations have the form 

(a\T" b (z)\b) = (a\V b \b)+J2 J (a\c) (c\V c R c (z)\c) {c\T cb (z)\b) (8.15) 

where (a|c) are Poincare group Racah coefficients, which are the unitary transformation that relate three-body 
Poincare irreducible bases constructed using pairwisc coupling in different orders. These coefficients, which have 
the form 

(a\c) = 6(p - p')S^S(m - m')6 n ,n m ^(d a , d c ), (8.16) 

with d a and db distinct sets of invariant degeneracy parameter, replace the non-relativistic permutation operators. 

The construction of the kernel is facilitated by the fact that the two-body eigenfunctions of (|8.4[) and (|8.5[) arc 
identical. The kernel of the relativistic Faddecv equation can be directly related to the non-relativistic two-body t 
using the following relations: 

(c'\V c R c (z c )\c) = (c'\T c (z c )(z c - Mo)" 1 ^) = 
(c'\V c \c-)(z c - Mo)- 1 = (c'\M c - M \c-)(z c - Mo)- 1 = 

; — , , p „ T, , ,, 9 =(c\t nr {k c )\c)(z c - M ) 8.17 

uJi = ^/fc 2 + m? , (8.18) 



where 



which holds for the half shell kernel; this can be used as input to construct the fully off-shell kernel using the first 
resolvent identity [5(J 

W) = T c (z c ) + T c (/) | 4gL T c fe) + z c (8.19) 

Alternatively, this kernel has also been computed using an iterative procedure based on a non-linear integral 



equation[51|. 



The differences with the non-relativistic three-body calculations are the different off-shell dependence dictated by 
(|8.17p . the differences in the Poincare group Racah coefficients (|8.16|) and the non-relativistic permutation operators. 
These differences show up for the first time in the three-body system, since our two-body interactions are designed 
to reproduce the same experimental two-body cross sections as the non-relativistic calculations. 

Solving these equations leads to three types of predictions: binding energies, 

M|*) = A|#) (8.20) 
\9) = E(V)\9u) MuY&bt) = A|* w ), (8.21) 



scattering probabilities (N — 3 only), 

' /in '/ \ *i /I ~ W^btfl^bti 



%| 2 H(*|l*r>| 2 H<^l* te )| 2 , (8-22) 



electromagnetic and weak current matrix elements 

(Vfinom) = (* f \Ai»(o)Atm , (8.23) 

where 1^(0) is a current that is conserved, covariant, and clusters in the representation f|6 . 1 T[) of the three-body 
dynamics. 

In what follows wc discuss three applications of Poincare invariant quantum mechanics that illustrate its ability to 
model a variety reactions where relativity may be important. 
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A. Relativisitic spin rotations in low energy A y 



A calculation by Miller and Schewnk [52| suggested that Wigner rotations might have an observable effect on the 
polarization observable A y for low-energy p-d scattering. Comparison of three-body calculations based on Poincarc 
invariant quantum mechanics [34| and non-relativistic quantum mechanics using the same realistic CD-Bonn inter- 
action [49| as input indeed show a surprising sensitivity of A y to Wigner rotations. These calculations, which arc 
shown in Fig. 8.1, compare the non-relativistic result (dotted curve), the relativistic result without Wigner rotations 
(dashed curve) and the full relativistic calculation (solid curve) to data [HI, HH . While the relativistic effects move 
the calculations away from the data, this calculation illustrates that the relativistic effects cannot be ignored in these 
calculations, even at these low energies. 



B. Relativistic Effects in Exclusive pd Breakup 



The value of Poincare invariant quantum mechanics is that it provides a consistent framework to study strong 
interactions in the few GeV energy scale. At this scale it is more efficient to perform calculations using direct 
integration [HBL [56j rather than with partial wave expansions. The feasibility of using Poinca'e invariant quantum 
mechanics to treat nucleon deuteron-scattering at these energy scales was established by solving the Faddeev equation 
of section ?? using Malflict-Tjon [57| interactions to model the nucleon-nucleon potential. The two-body interactions 
were included using the method discussed above. Convergence of the solutions of the Faddeev equations was established 
up to 2GeV |3ll - l33l |. In three-body reactions there are many observables that can be used to test the sensitivity of 
relativistic effects. One interesting observable is the cross section when the outgoing protons in a breakup reaction 
are measured at symmetric angles relative to the beam direction. These cross sections were computed |32| in non- 
relativistic and Poincare invariant three-body models using the same Malfliet-Tjon two-body interactions as input. 

Fig. 8.3 shows cross sections for different choices of angles symmetric about the beam direction. The solid curve is 
the relativistic calculation while the long dashed curve is the non-relativistic one. The other two curves compare the 
exact calculation to the first terms in the multiple scattering series both for the relativistic and non-rclativistic cases. 
As the angle is increased the relativistic and non-rclativistic curves, exhibit different behavior. For this kinematic 
configuration the multiple scattering series converges quickly, although this result depends on what is measured. 
Fig. 8.2 shows similar plots for non-symmetric angles. Again the first order multiple scattering calculations work 
reasonably well and there is a definite difference between the relativistic and non-relativistic predictions. In both 
cases the data [HI, has the same qualitative behavior as the relativistic calculations. 



C. Exchange currents in electron-deuteron scattering 

The last application involves electron scattering off of nuclear targets at values of momentum transfer Q 2 appropriate 
to J-lab experiments. In Poincare invariant quantum mechanics electron scattering observables in the one-photon- 
exchange approximation can be expressed in terms of matrix elements of a conserved covariant current /^(x) which 
should have a cluster expansion 

to = E to + E TO + E W + ■■■■ (s- 24 ) 

Both Poincare covariance, current conservation, and cluster properties put dynamical constraints on the current 
operator. 

The deuteron is the simplest electromgnetic target that is sensitive to the two-body part of the current. While 
a general method for constructing I^{x) based on dynamical considerations is not known, the constraints can be 
satisfied by using the Wigncr-Eckart theorem for the Poincarc group, which amounts to computing a maximal set of 
linearly independent current matrix elements and using covariance and current conservation to generate the remaining 
matrix elements. Different model two-body currents can be tested in this framework. For clastic scattering off of 
a deuteron there are three independent observables which can be taken as, A(Q 2 ), B(Q 2 ), and T2q(Q 2 ,70°). The 
input to a calculation is a deuteron wave function, a dynamical representation of the Poincare group, nucleon form 
factors [59M63T ]. and a model exchange current [(34|. The calculations illustrated in Figs. 8.4-8.6 use a model of the 
deuteron with a light front kinematic symmetry. The dynamical representation of the Poincare group is constructed 
from the Argonne V18 inter action [48j], and the exchange current is the long-range part of a "pair current" derived 
from the one-pion-exchange part of the V18 interaction. Figs. 8.4, 8.5, and 8.6 show comparisons of A(Q 2 ), B(Q 2 ), 
and T2o(Q 2 , 70°) to experimental data with and without the exchange current. Two different implementations of the 
Poincare group Wigncr-Eckart theorem are responsible for the small difference in the curves labeled II and III. 
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FIG. 1: Comparison of relativistic and non-relativistic calculations of the observable A y at low energies. 



These three calculations illustrate both the power and flexibility of Poincae invariant quantum mechanics as a tool 
to study systems of str ongly i nteracting particles at scales up to a few GcV. Data shown for A are from (65l - [75j . for 
B are from [H |zj [zl, [2H73 , and for T 20 from (ziil- 

These calculations demonstrate that Poincare invariant quantum mechanics is a useful framework for making 
realistic models of system of strongly interacting particles at the few-GeV energy scale. Some of these effects extend 
to surprisingly low energies. 

This work was performed under the auspices of the U. S. Department of Energy, Office of Nuclear Physics, under 
contract No. DE-FG02-86ER40286. 
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FIG. 2: Comparison of rclativistic and non-relativistic calculations of exclusive proton deuteron breakup scattering at non- 
symmetric angles. 
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FIG. 3: Comparison of relativistic and non-relativistic calculations of exclusive proton deuteron breakup scattering at non- 
symmetric angles. 





T 9n (Q 2 ,70°): BBBA; II, III 




FIG. 6: Comparison of a front-form calculation of the elastic electron-deuteron tensor polarization with and without "pair 
current" contributions. 
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